Automatic implementation of material laws: Jacobian calculation in a finite element code with TAPENADE
نویسندگان
چکیده
In an effort to increase the versatility of finite element codes, we explore the possibility of automatically creating the Jacobian matrix necessary for the gradient-based solution of nonlinear systems of equations. Particularly, we aim to assess the feasibility of employing the automatic differentiation tool TAPENADE for this purpose on a large Fortran codebase that is the result of many years of continuous development. As a starting point we will describe the special structure of finite element codes and the implications that this code design carries for an efficient calculation of the Jacobian matrix. We will also propose a first approach towards improving the efficiency of such a method. Finally, we will present a functioning method for the automatic implementation of the Jacobian calculation in a finite element software, but will also point out important shortcomings that will have to be addressed in the future.
منابع مشابه
Fast Automatic Differentiation Jacobians by Compact LU Factorization
For a vector function coded without branches or loops, a code for the Jacobian is generated by interpreting Griewank and Reese’s vertex elimination as Gaussian elimination and implementing this as compact LU factorization. Tests on several platforms show such a code is typically 4 to 20 times faster than that produced by tools such as Adifor, Tamc, or Tapenade, on average significantly faster t...
متن کاملFast Ad Jacobians by Compact Lu Factorization
For a vector function, coded without branches or loops, code for the Jacobian is generated by interpreting Griewank and Reese’s Vertex Elimination as Gaussian elimination, and implementing this as compact LU factorization. Tests on several platforms show such code is typically 4 to 20 times faster than that produced by tools such as Adifor, Tamc or Tapenade, on average significantly faster than...
متن کاملBuckling Analysis of Rectangular Functionally Graded Plates with an Elliptic Hole Under Thermal Loads
This paper presents thermal buckling analysis of rectangular functionally graded plates (FG plates) with an eccentrically located elliptic cutout. The plate governing equations derived by the first order shear deformation theory (FSDT) and finite element formulation is developed to analyze the plate behavior subjected to a uniform temperature rise across plate thickness. It is assumed that the ...
متن کاملFINGER: A Symbolic System for Automatic Generation of Numerical Programs in Finite Element Analysis
FINGER iS a Lisp-based system to derive formulas needed in finite element analysis, and to generate FORTRAN code from these formulas. The generated programs can be used with existing, FORTRAN-based finite element analysis packages. This approach aims to replace tedious hand computations that are time consuming and error prone. The design and implementation of FINGER are presented. Techniques fo...
متن کاملForward-Mode Automatic Differentiation in Julia
We present ForwardDiff, a Julia package for forward-mode automatic differentiation (AD) featuring performance competitive with low-level languages like C++. Unlike recently developed AD tools in other popular high-level languages such as Python and MATLAB, ForwardDiff takes advantage of just-in-time (JIT) compilation to transparently recompile AD-unaware user code, enabling efficient support fo...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Computers & Mathematics with Applications
دوره 72 شماره
صفحات -
تاریخ انتشار 2016